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Abstract 

Objective: This study presents and evaluates a semi-automated algorithm for quantifying 
malalignment in complex femoral shaft fractures from a single intraoperative cone-beam CT 
(CBCT) image of the fractured limb. 

Methods: CBCT images were acquired of complex comminuted diaphyseal fractures created in 
9 cadaveric femora (27 cases). Scans were segmented using intensity-based thresholding, 
yielding image stacks of the proximal, distal and comminuted bone. Semi-deformable and rigid 
affine registrations to an intact femur atlas (synthetic or cadaveric-based) were performed to 
transform the distal fragment to its neutral alignment. Leg length was calculated from the 
volume of bone within the comminution fragment. The transformations were compared to the 
physical input malalignments. 

Results: Using the synthetic atlas, translations were within 1.71 ±1.08 mm (medial/lateral) and 
2.24 ±2.11 mm (anterior/posterior). The varus/valgus, flexion/extension and periaxial rotation 
errors were 3.45 ± 2.6°, 1 .86 ± 1 .5° and 3.4 ± 2.0°, respectively. The cadaveric-based atlas yielded 
similar results in medial/lateral and anterior/posterior translation (1.73 ±1.28 mm and 
2.1 5 ±2.1 3 mm, respectively). Varus/valgus, flexion/extension and periaxial rotation errors 
were 2.3 ±1.3°, 2.0 ±1.6° and 3.4 ±2.0°, respectively. Leg length errors were 1.41 ±1.01 mm 
(synthetic) and 1.26 ±0.94 mm (cadaveric). The cadaveric model demonstrated a small 
improvement In flexion/extension and the synthetic atlas performed slightly faster 
(6 min 24 s ± 50 s versus 8 min 42 s ± 2 min 25 s). 

Conclusions: This atlas-based algorithm quantified malalignment in complex femoral shaft 
fractures within clinical tolerances from a single CBCT image of the fractured limb. 
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Introduction 

Femoral shaft fractures are commonly encountered in ortho- 
paedic trauma and occur most frequently as a result of high- 
energy trauma in individuals under the age of 35 [1]. Closed 
intramedullary (IM) nailing of the femur has become the 
standard of care for adults due to its high union rates, 
low complication rates and minimally invasive nature [2]. 
Clinically, in complex comminuted fracture reduction, the two 
main fragments are normally aligned in anterior-posterior 
and medial-lateral translation and the angle between the 
long axes of the bones is then restored. Leg length and 
periaxial rotation are restored about a guidewire and nail. 
Fluoroscopic guidance enables the visualization of the 
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fracture fragments; however, the resolution of these images 
and their 2D nature limit the ability to successfully align 
complex fractures in all six degrees of freedom (DOF). 

Despite its popularity, complications relating to malalign- 
ment occur in IM nailing [3-5]. Current clinical tolerances for 
alignment of femoral shaft fractures are 5-10° varus/valgus 
and 10° flexion/extension angulation, 15° periaxial rotation, 
and 1.0-1.5 cm leg length and medial-lateral and anterior- 
posterior translation [6-9]. In some series, severe malalign- 
ments in excess of these tolerances have been reported in 
19-30% of patients, leading to malunions [2, 3, 10-12]. 
Fluoroscopy-based workflows combined with navigation have 
been used in femoral shaft fracture reduction with IM nailing 
[13, 14]. Malalignment is determined with respect to 
landmarks on the contralateral limb and through manual 
digitization of the fracture fragments. These systems have 
achieved only variable results with respect to periaxial 
rotation [13, 15]. Three-dimensional imaging in the form of 
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cone-beam CT (CBCT) paired with flat-panel detectors can 
provide volumetric data from a series of projections obtained 
in a single rotation of a unit about a patient [16]. CBCT offers 
the ability to obtain 3D data with sub-millimetric isotropic 
resolution and a large field of view. CBCT images enable the 
use of surface features and 3D-based analyses to guide 
fracture alignment. 

The objective of this work was to develop and validate 
a semi-automated algorithm in which malalignment 
between superior and inferior fracture fragments is deter- 
mined based on a single CBCT image of the injured 
femur. We hypothesized that an atlas-based approach would 
enable the quantification of the malalignment within clinical 
tolerances in 6 DOF, even in the presence of severe 
comminution. This study further aimed to compare the 
performance of a synthetic atlas (based on averaged anatom- 
ical data) with that of an atlas generated based on a single 
human femur. 

Materials and methods 

Atlas construction 

Two atlases were created from CT scans performed on a 
GE LightSpeed VCT scanner (GE Corp., Fairfield, CT): 
a right intact synthetic femur (Pacific Research Laboratories 
Inc., Vashon, WA) and a left intact cadaveric femur 
(from a 50-year-old male with no visible pathology). 
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both at a voxel size of 0.67 x 0.67 x 1.25 mm. 
Each atlas scan was imported into an image-analysis 
software platform (AmiraDEV 5.2, Visage Imaging, Berlin, 
Germany), segmented using intensity-based thresholding, and 
saved as an image stack that included only the voxels 
containing bone. 

Each atlas was aligned so that the main anatomical axes 
were coincident with the Cartesian coordinates of the global 
image analysis system. A module within the AmiraDEV 
environment was used to calculate the three main axes of 
mass distribution; the axis with the smallest moment of inertia 
was defined as the z-axis of the atlas shaft and aligned to the 
global z-axis [17]. Landmarks were then placed on the most 
posterior point of the femoral condyles using coronal slices 
and this axis was aligned with the global x-axis, again using 
the principal axes (Figure 1). The medial-lateral (x), anterior- 
posterior (y) and longitudinal (z) translation, the flexion/ 
extension (a) and varus/valgus (P) angulation, and the 
periaxial ((/?) rotation angles could thus be used to describe 
malalignment in clinically relevant combinations of transla- 
tions and rotations. 

Fracture model and scan acquisition 

Nine dried cadaveric left femora with no visible pathology 
(Pacific Research Laboratories Inc.; no demographic data was 
available from the supplier) were obtained, stripped of soft 
tissue and mounted in a mechanical jig (Figure 2). The jig was 



Figure 1. Isosurfaces of the synthetic (left) 
and cadaveric (right) atlases with the 6 DOF 
of malaUgnment labelled. Surface features on 
the synthetic atlas appear much smoother and 
less pronounced compared to those of the 
cadaveric atlas. 
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composed of a clamp mounted on a ball joint that enabled 
movement in 6 DOF to align the sample. The specimen 
was then secured at its distal end to a platform that had 
been mounted on a series of six mechanical stages 
(Newport Corporation, Irvine, CA) permitting translation 
(accurate to 0.05 mm) and rotation (accurate to 0.5°) about the 
X-, y- and z-axes. 

Within the jig, a diaphyseal osteotomy was performed on 
each specimen, removing bone of various lengths (from 7.5 to 
76 mm) to simulate fracture comminution. Removed segments 
of bone were measured using calipers accurate to 0.05 mm. 
Examination of markings placed on the bone ensured the two 
fragments initially maintained their neutral alignment pos- 
itions within the jig. The comminuted bone fragment was 
malpositioned in the fracture gap. The jig was placed within 
the field of a CBCT scanner (CB MercuRay scanner, Hitachi 
Medical Systems America, Twinsburg, OH) and CBCT 
images were acquired of the specimens under random 
malalignments (voxel size 0.377 x 0.377 x 0.377 mm; field 
of view 30.5 cm; 80kVp and 10 mA). Each bone was placed 
in three different configurations, leading to a total of 
27 malaligned scans of 9 femora with 9 different fracture 
gaps (ranges are summarized in Table I). Longitudinal 
translation {z) of the distal fragment was performed to 
vary the length of the proximal and distal fragments that 
were visible within the field of view. Images were seg- 
mented using intensity-based thresholding and saved into 
three image stacks: the proximal/superior fragment, the 



Table I. Ranges of medial-lateral translation {x), anterior-posterior 
translation (y), flexion/extension angulation (a), varus/valgus angulation 
(/3), periaxial rotation (0), and leg length (z) used to create the 
malalignment conditions. 



Parameter 


Min. value 


Max. value 


X (mm) 


0.80 


68.55 


)' (mm) 


0.35 


54.30 


z (ram) 


7.5 


76.0 


an 


0 


13.5 


/3(°) 


0 


13.5 




3.5 


83.0 



comminution or fractured volume, and the distal/inferior 
fragment. 

Semi-automated registration 

A semi-automated algorithm was written in Tool Command 
Language for AmiraDEV (run on an Intel Xenon 2 CPU 
2.33 GHz with 8 GB RAM). Figure 3 is a flow chart outlining 
the main steps within the registration algorithm. The three 
image stacks were loaded into the software and the user 
was prompted to specify the fracture side. If the fractured 
femur was opposite to that of the atlas (left versus right), the 
atlas was mirrored about the z-axis. Isosurfaces of 
the fractured femur and the atlas were visualized in 3D and 
the atlas was roughly spatially aligned to the proximal 
fragment using a manual, interactive transform editor 



DOI: 10.3109/10929088.2014.894126 



CT algorithm for femoral fracture reduction 51 



op S 




■3 -S .c 



5 
o 

E 



m .X 



-J 



52 M. C. Crookshank et al. 



Comput Aided Surg, 2014; 19(1-3): 48-56 



(9 DOF, allowing scaling, translational and rotational adjust- 
ments) to provide an initial estimate of its position with 
respect to the proximal fragment. Shearing was not included 
in the transformation of the atlas to ensure that the relation- 
ship between the principal axes and the global coordinate 
system remained constant. An interactive editor was used to 
crop the atlas so that it contained only the region overlapping 
the proximal fragment. The cropping step was included so 
that the affine registration did not attempt to align anatomical 
landmarks to areas in the full atlas scan not present in the scan 
of the fractured specimen and to decrease the time required 
for the registration. 

Following cropping of the atlas, an automated affine 
registration with 9 DOF was performed allowing translation, 
rotation and anisoscaling of the atlas to best align it to the 
superior fragment [17]. The registration initially down- 
sampled both image sets by 8x in all directions to decrease 
the computational time required for the coarser transforms. 
Both images were then sequentially upsampled to perform 
finer transform steps until their initial resolution was 
restored, with a step length of 0.01 mm for fine tuning. The 
registration used a quasi-Newtonian optimization method to 
determine the final registration, relying on a normalized 
mutual information metric which was used to account for the 
differences in the image intensities arising from different 
imaging modalities, CBCT and CT, and resolutions [18, 19]. 
Only voxels that had been included in the initial segmentation 
of both the atlas and the superior fragment were included 
in the registration. 

The inverse of the registration transform without the 
scaling factors was then calculated and applied to all image 
stacks (the full atlas, the proximal and distal fracture 
fragments, and the comminution volume). This aligned the 



Figure 4. (a) The synthetic atlas, proximal 
and distal fracture fragments in their scanned 
positions, (b) The atlas, registered to the 
proximal fracture fragment, with the inferior 
fragment in its malaligned position, (c) The 
proximal fragment, with the registered atlas, 
and the aligned distal fracture fragment. 




atlas and proximal fragment with the global coordinate 
system. The center of the bone contained within the most 
distal slice (user-selected) of the proximal fragment was 
calculated based on its voxel intensities and its distance to the 
global origin was determined. The proximal and distal 
fragments, atlas and comminution volumes were translated 
so that the center of the most distal slice was positioned at 
the origin. As such, the global origin, about which the 
transformation for the registration of the distal fragment 
occurs, would be identical to the most distal point of the 
proximal fragment, about which such a translation would be 
performed clinically. The distal fracture fragment was 
then registered to the atlas using rough manual alignment 
followed by a rigid affine registration with 6 DOF (allowing 
translation and rotation of the distal fragment but no scaling) 
(Figure 4). 

The algorithm next prompted the user to select the slices of 
the proximal and distal fracture fragments that bordered the 
fracture gap but contained full cross-sections of bone 
(Figure 5). The cross-sectional area of bone (bone-containing 
voxels) was averaged for the proximal and distal slices and 
used to determine the longitudinal translation (z) by dividing 
this area by the total amount of bone in the comminuted 
volume. 

The transformation matrix from the registration of the 
distal fragment to the atlas was then decomposed into the 
corresponding Euler angles and translations. The order of 
decomposition was based on the typical clinical approach and, 
therefore, was performed as rotations about the x-, y- and, 
finally, z-axes. This was performed using Equation Set 1, 
where T is the 4 x 4 3D transformation matrix yielded from 
the affine registration [20]. Due to the initial alignment of the 
atlas to the global origin, this corresponded to the a, P and (f> 




DOI: 10.3109/10929088.2014.894126 



CT algorithm for femoral fracture reduction 53 




Table II. The average difference between the repeated measures for 
medial/lateral (M/L) translation, anterior/posterior (A/P) translation, 
leg length, flexion/extension (F/E) angulation, varus/valgus (V/V) 
angulation and periaxial rotation. The maximum difference between 
the repeated measurements is also indicated, along with the fracture gap 
at which the maximum difference occurred. 




Figure 5. A schematic indicating the cross-sectional areas of bone (black 
and white hatching) in the gap-bordering slices of the proximal (Al) and 
distal (A2) fragments. The comminution fragment (V) is also shown to 
the left of the atlas. 



that were displayed for the user, along with the volume-based 
leg length correction. 
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Performance evaluation 



The calculated values were compared to the input values from 
the precision stages. Differences in the two values were 
calculated for each scan using both atlases, and the error in 
each of the six alignment measurements was compared 
between the two atlas models using a two-tailed Student's 
paired t test with a significance level of 0.05. The time 
required for the algorithm to run was also recorded and 
compared between the two atlas models using a two-tailed 
Student's paired t test with a significance level of 0.05. To 
determine whether the volume of missing bone and the 
amount of bone visible in the scan effected the accuracy of the 
algorithm, a Pearson correlation was performed between each 
of the six alignment measurements and the gap length and 
minimal amount of bone visible in either the proximal or 



Synthetic atlas 



Cadaveric atlas 



Parameter 



Mean Max. Mean Max. 

error error (gap size) error error (gap size) 



M/L translation (mm) 0.81 2.24 (60.2 ram) 1.52 2.27 (60.2 mm) 



A/P translation (mm) 1.36 

Leg length (mm) 0.25 

F/E angulation (°) 1.2 

V/V angulation (°) 0.9 

Periaxial rotation (°) 1.8 



3.63 (62.5 ram) 
0.44 (9.9 mm) 
3.3 (28.6 ram) 
2.3 (62.5 ram) 
3.2 (28.6 ram) 



1.54 

0.36 

2.0 

1.2 

1.4 



5.97 (62.5 mm) 
1.12 (62.5 mm) 
3.3 (76.0 mra) 

2.5 (39.1 mra) 

2.6 (39.1 mra) 



distal fragment in the scan (the length of the shortest fracture 
segment). Finally, to determine repeatability, the semi- 
automated algorithm was run a second time for both atlas 
models on one scan from each bone (9 scans in total with 
different gap lengths) and the absolute difference between the 
first and second measurements was calculated. 

Results 

The algorithm was able to calculate the relative malalignment 
within current clinical tolerances in all 6 DOF using both the 
synthetic and cadaver-based atlas models. Table II summar- 
izes the results from the 9 repeated measurements for each of 
the 6 DOF. The algorithm was found to be highly repeatable, 
with an average difference between the initial measurement 
and the repeated measurement of < 1.55 mm for all transla- 
tions and 2° for rotations and angulations. The paired 
Student's t test showed no significant difference for the 
accuracy of the quantification of the A/P and M/L translation, 
leg length, varus/valgus angulation and periaxial rotation 
(/7>0.36 for all of these). However, the flexion/extension 
angulation calculated based on the cadaver-based atlas was 
found to be more accurate than that calculated based on the 
synthetic femur, with a mean difference of 1.2° between the 
two methods [p = 0.043). 

The average time required for the algorithm to calculate 
the malalignment was 6:24 min (±50 s, maximum 8:15min) 
and 8:42 min (±2:25 min, maximum 13:48 min) for the 
synthetic-based and cadaver-based atlases, respectively. 
From the paired Student's t test, the cadaver-based atlas 
model required significantly more time to perform the 
registration and calculate the malalignment, with a mean 
difference of 2: 18 min (/?< 0.001). The amount of time 
required by the user did not differ between the two atlases and 
was approximately 3 minutes. 

The Pearson correlation indicated that the size of the 
fracture gap had a positive correlation with the error in the 
calculation of the M/L and A/P translation (r = 0.42, 
slope = 0.019 and r=0.66, slope = 0.058, respectively) as 
well as the calculated errors in flexion/extension angulation 
(r = 0.55, slope = 0.058) for the synthetic -based atlas and the 
A/P translation and periaxial rotation (r = 0.68, slope = 0.060 
and r=0.49, slope = 0.042, respectively) for the cadaver- 
based atlas (all significant at the 0.05 level). No significant 
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correlation was found between the other parameters for either 
atlas and the size of the fracture gap. The length of the 
shortest fracture segment was found to have a negative 
correlation with the error in the calculation of the flexion/ 
extension angulation and periaxial rotation for the synthetic- 
based atlas (r=:— 0.56, slope—— 0.10 and r=— 0.50, slope 
= —0.069, respectively) at a level of significance of 0.01. 
No significant correlation was found between the length of the 
shortest fracture segment and the accuracy of the other 
parameters in the synthetic-based atlas or the cadaver-based 
atlas. 

Discussion 

This study presents a semi-automated algorithm that can 
calculate the relative alignment of the two main fragments in 
femoral shaft fractures with comminution in 6 DOF within the 
current clinical tolerances [6-9, 21-23]. Two distinct atlases 
were investigated for use in the algorithm: a synthetic-based 
atlas (representing a population average) and a cadaveric- 
based atlas generated from a single specimen. In both cases, 
the algorithm was able to quantify the degree of malalignment 
based on a single CBCT image of the fractured femoral shaft. 
The algorithm was performed in two main steps: a rigid 
registration with anisoscaling of the atlas to the proximal 
fracture fragment, followed by the rigid registration of the 
distal segment to the atlas, yielding the relative alignment 
between the two segments. 

The accuracy of some of the alignment parameters was 
affected by the amount of bone that was contained in the 
comminution volume and the amount of bone that was visible 
within the scan (the fragment length). This was consistent 
with our expectation, as both shorter bone fragments and 
increased comminution length lead to less information being 
available with which the algorithm can calculate the align- 
ment. However, the size of this effect was small, with a 
maximum error of 0.06 mm and 0.06° per mm of additional 
comminution and 0.10° per mm of the minimum visible 
fragment length. Also, the accuracy of the majority of 
parameters was found to have no significant correlation to the 
size of the fracture gap or minimum fragment length. For 
fracture gaps of up to 76 mm, and with a minimum visible 
fragment length of 42 mm, the algorithm was able to quantify 
successfully malalignment of the fracture in 6 DOF to current 
clinical standards. 

The cadaveric-based atlas yielded slightly more accurate 
results (1.2°) in the calculation of the flexion/extension 
angulation, with no significant differences being found 
between the two atlases in the other 5 DOF. The cadaveric- 
based atlas required significantly more time to perform the 
registration and alignment quantification, probably due to the 
presence of more varied features along its surface. As the 
synthetic-based atlas represents a population average, some of 
the surface features, such as the linea aspera, are less defined 
and there are no irregular surface features. Additionally, the 
synthetic-based atlas has a more uniform density, yielding 
more uniform segmentation and a smoother surface. Future 
work may investigate the development of an optimized atlas 
incorporating features of the synthetic and cadaveric femora 
to provide improved accuracy and computational speed. 



This algorithm was not fully automated in that it requires 
initial manual alignment of the atlas to the proximal fragment 
and of the distal fragment to the atlas. The manual input is 
simple and rapid, requiring approximately 3 minutes in total. 
The robust and repeatable results from the algorithm are 
encouraging when considering its potential for translation to 
the clinical setting. Current navigation systems require 
acquisition of multiple fluoroscopic images on both the 
fractured and intact limb and the selection of anatomical 
landmarks, from which the quantification of leg length and 
periaxial rotation are calculated. Although time requirements 
for the use of fluoroscopy-based navigation have not been 
widely documented, Kendoff et al. reported times from 
16-18 min [13]. Other navigated methods have been reported 
to increase overall procedural time, despite a reduction in 
radiation exposure time, due to the image acquisition and 
registration process, as well as the additional set-up and 
equipment required. The present algorithm could be incorpo- 
rated into current tracking systems using a single CBCT scan 
of the fractured limb (with mounted distal and proximal 
trackers). 

The presented algorithm also provides quantification of the 
fracture in 6 DOF, which is not present in some of the current 
navigation techniques which solely provide quantification of 
the leg length and periaxial malrotation, and rely on surgeon 
judgment for the additional translational and rotational 
alignments [14, 24, 25]. Additionally, the single intraoperative 
CBCT scan could reduce radiation exposure for surgical staff 
(who could be shielded during image acquisition) with no 
additional fluoroscopic imaging being needed for reduction. 

Current fluoroscopy-based navigation techniques have 
several inherent limitations. Many require images of an 
intact contralateral limb (which is not always available). 
Additionally, the metrics used to evaluate the current methods 
are 2D-based and, therefore, may have systemic errors leading 
to the varied reports of accuracy, particularly in the 
determination of periaxial rotation. Hofstetter et al. [14], in 
one of the initial reports on using fluoroscopic landmark 
identification for the determination of periaxial rotation, 
found an average accuracy of 3.6° with a maximum of 9° 
difference in unfractured pairs of femora. Weil et al. [26] 
reported average accuracies of 1.9 mm for leg length calcu- 
lations and 2.5°, using a pre-osteotomy scan of the same 
femur for a fracture gap of 50 mm. However, these results 
were based on a single femur and, instead of using the 
contralateral limb for length and periaxial rotation determin- 
ation, images of the same femur from before the fracture was 
simulated were used during the matching scheme. In a series 
of three patients, Kendoff et al. [13] measured anteversion 
postoperatively using CT scans and found, in three cases, that 
the periaxial rotation was restored to within 4° of that of the 
healthy side. Due to the reliance on the presence of an intact 
contralateral limb, such fluoroscopic methods could not be 
used to determine periaxial rotation or leg length in patients 
with bilateral fractures. Furthermore, no information was 
provided on the level of comminution in the three treated 
fractures, making evaluation of the alignment difficult to 
assess. Keast-Butler et al. [25] observed no significant 
difference between the reduction accuracy achieved with 
fluoroscopy-based navigation and conventional fluoroscopy 
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(maximum error was 20° for both groups). Citak et al. [15] 
found that, when using non-orthogonal images, with an angle 
of less than 50° between the two images, the periaxial rotation 
was measured to be from 2.3 to 7.7° different from that 
calculated from a CT-based measurement of anteversion of 
the same femur. Therefore, off-angle fluoroscopic images can 
further decrease the calculation of the periaxial rotation. 

Previous computer assisted methods based on CT and 
fluoroscopic guidance have mainly focused on determination 
of periaxial rotation [13, 15, 26-32]. The FRACAS system 
developed by Joskowicz et al. [30] is CT-based and enables 
3D visualization of the fracture fragments, but requires 
preoperative CT imaging of the fractured and contralateral 
healthy limb and their registration to fluoroscopic images. 
In pilot studies [31], landmarks placed on the posterior 
femoral condyles and the neck axis of the healthy and 
fractured side were able to guide the fracture to within 1.0 to 
4.5° of the alignment of the healthy limb, based on the 
measurement of the selected landmarks. However, this 
method requires substantial radiation exposure through pre- 
operative CT and intraoperative fluoroscopy. 

Three-dimensional CT-based surface information has been 
used, by our group and others, to determine periaxial rotation 
of fractured femurs with accuracies from 1.6 to 3° [33, 34]. 
Each group used a surface-matching approach and resampled 
the femoral shaft about its center to yield surface maps of the 
femur. Additionally, Moghari and Abolmaesumi [35] obtained 
promising results using a statistical atlas for the registration of 
bone fragments in femoral neck fractures, reporting registra- 
tion errors of approximately 1.6 mm. These approaches are 
advantageous as they require no information from the 
contralateral limb. A surface matching approach could 
ultimately be incorporated into the atlas-based approach to 
provide fine tuning of the periaxial rotation. 

More recent work has integrated surface matching with a 
principal axis-based approach to allow the determination of 
malalignment of fractured femora in 6 DOF [36]. This work 
demonstrated robust performance that can be directly 
compared to the current algorithm (the principal axis/ 
surface-based approach was applied to the same data set of 
CBCT images). The principal axis/surface-based approach 
achieved clinical tolerances in all but two cases (one case at 
the extreme of comminution length and one case with the 
minimum length of fracture fragment), while the present work 
achieved clinical tolerances in all 27 fracture cases. A paired 
Student's t test showed no significant difference between the 
principal axis/surface-approach and the atlas-based method in 
the present work for all alignment parameters {p > 0.20) but 
periaxial rotation (/?< 0.001; both atlases achieved greater 
accuracy) and the flexion/extension angulation (p — 0.02; the 
cadaveric-based atlas achieved greater accuracy). However, 
the principal axis/surface method was more rapid (average 
time for analysis of 5min58s) than the atlas-based method 
presented here (average analysis times of 6min24s and 
8min42s for the synthetic and cadaver-based atlases, 
respectively). 

Intensity-based threshold segmentation of cortical bone is 
robust even in the presence of soft tissue, particularly at the 
high resolution and image quality that can be achieved with 
CBCT. As such, the lack of soft tissue in the CBCT images in 
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this study would likely have little effect on the segmentation 
or algorithm performance. Additionally, the algorithm cur- 
rently requires the input of separate image stacks for the 
proximal, distal and comminution bone volumes. The high 
quality and contrast in the CBCT images, the femoral shaft 
geometry and the absence of trabecular bone make the 
intensity-based segmentation and cropping of the relevant 
bony volumes straightforward and repeatable. 

The algorithm in this study was not tested to its limits for 
fracture gap, minimal length of bone or fracture malalign- 
ment. However, alignments in excess of those that were 
created using the precision stages in this study would be less 
relevant to clinical practice, as they would be readily apparent 
on 2D images and may be partially reduced by the placement 
of a patient in traction prior to surgery. The fracture gap of 
76 mm represents substantial bone loss. Fracture gaps in 
excess of this length would likely start to approach the 
epiphyses and, additionally, would limit the amount of bone 
visible within the field of view. Future analyses will 
determine the minimum fracture fragment length necessary 
for accurate registration and alignment. 

Overall, this algorithm enables repeatable and robust 
quantification of malalignment of femoral shaft fractures in 
6 DOF within clinical tolerances based on a single initial 
intraoperative CBCT image of the fractured limb. The 
incorporation of this algorithm with a navigation system 
that can support intraoperative CBCT images presents a 
promising alternative to the current fluoroscopy-based navi- 
gation. Integration and validation of 3D algorithms in 
intraoperative CBCT imaging/navigation systems and their 
associated workflow will ultimately demonstrate their trans- 
lational potential. 
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